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On  Thermodynamics  and  the  Nature  of  the  Second  Law 
for  Mixtures  of  Interacting  Continue 


by 


A.  E.  Green"''  and  P.  M.  Naghdi* 


Abstract . This  paper,  which  may  be  regarded  as  a continuation  of  a previous 
paper  on  thermodynamics  of  single  phase  continue,  is  concerned  with  a new 
approach  to  thermomechanics  of  multiphase  continue  and  extends  the  previous 
ideas  and  procedure  (Green  and  Naghdi  1977)  to  mixtures  of  interacting  continua. 

In  particular,  it  contains  (a)  a proposal  of  a new  approach  for  obtaining 
restrictions  on  constitutive  equations,  (b)  an  appropriate  mathematical  state- 
ment of  the  second  law  for  mixtures  and  (c)  the  nature  of  restrictions  placed 
by  the  latter  on  constitutive  results  representing  the  thermo-mechanical 
behavior  cx  mixtures  with  different  constituent  temperatures.  Our  point  of 
departure  is  the  introduction  of  balances  of  entropy  and  the  use  of  a single 
energy  equation  for  the  whole  mixture  as  an  identity  for  all  motions  and  all 
temperature  distributions  after  the  elimination  of  the  external  fields.  This 
procedure  is  in  contrast  to  the  existing  approach  in  most  of  the  current  litera- 
ture on  continuum  theories  of  mixtures  based  on  the  use  of  a Clausius -Duhem 
type  inequality  (or  similar  entropy  inequalities)  for  mixtures.  Our  interpreta- 
tion of  the  second  law  is  similar  to  that  of  the  previous  paper  and  leads  us  to 
postulate  an  inequality  which  reflects  the  fact  that  for  every  process  associated 
with  a dissipative  mixture,  a part  of  the  external  mechanical  work  is  always 
converted  into  heat  and  this  cannot  be  withdrawn  from  the  mixture  as  mechanical 
work.  The  restriction  on  the  heat  conduction  vectors  is  considered  separately 
and  is  confined  to  equilibrium  cases  in  which  heat  flow  is  steady.  Also,  a 
restriction  on  the  specific  internal  energies  is  derived  when  the  mixture  is  in 
the  state  of  mechanical  equilibrium  with  all  its  constituents  at  a common 
spatially  homogeneous  temperature.  The  remainder  of  the  paper  deals  with  the 
constitutive  results  and  thermodynamic  restrictions  for  inviscid  fluids,  a 
fairly  detailed  consideration  of  the  problem  of  an  incompressible  Newtonian  viscous 
fluid  flowing  through  a rigid  porous  solid  whose  temperature  may  be  different 
from  that  of  the  fluid,  as  well  as  some  additional  remarks  on  the  implication 
of  the  use  of  an  entropy  inequality  of  the  Claus ius -Duhem  type  for  mixtures  as 
contrasted  with  the  thermodynamic  restrictions  resulting  from  the  procedure 
proposed  here  and  from  our  present  interpretation  of  the  second  law. 
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Introduction 


The  continuum  theory  of  mixtures  has  been  a subject  of  intensive  study  in 
recent  years  from  different  points  of  view;  and,  at  present, there  appears  to 
be  agreement  on  the  basic  thermodynamical  equations  of  the  theory.  Differences, 
however,  arise  over  the  manner  in  which  a second  law  of  thermodynamics  is 
interpreted  for  mixtures.  A second  law  is  usually  employed  to  place  restric- 
tions on  constitutive  equations;  and,  even  for  a single  phase  continuum,  this 
aspect  of  the  subject  is  a matter  of  some  controversy.  We  do  not  attempt  a 
detailed  review  of  the  literature  on  mixtures,  but  restrict  our  attention  to 
references  which  are  relevant  to  our  present  development.  Briefly,  an  account 
of  the  subject  with  extensive  references  prior  to  i960  is  given  by  Truesdell 
and  Toupin  (10).  An  alternative  approach  for  mixtures  with  a single  tempera- 
ture by  Green  and  Naghdi  (3 >5)  is  based  on  different  primitive  concepts  and 
consequently  some  of  their  basic  field  equations  are  of  different  forms  from 
those  of  Truesdell  and  Toupin  (10).  However,  it  has  been  shown  by  Green  and 
Naghdi  (5)  that  the  basic  equations  in  the  two  formulations  of  the  theory  are 
equivalent,  although  this  is  not  apparent  at  first  sight  in  the  case  of  some 
of  the  equations.  For  further  details  end  remarks  concernir^;  the  relation- 
ship between  the  two  forms  of  the  theory  and  the  interpretations  associated  with 
some  of  the  quantities  which  occur  in  the  equations,  we  refer  the  reader  to  Green 
and  Naghdi  (5,8)  and  to  a recent  review  articles  by  Atkin  and  Craine  (l). 

The  present  paper  is  a companion  to  the  previous  work  on  thermodynamics  of 
single  phase  continua  by  Green  and  Naghdi  (9a, b).  The  latter  paper  contains  a new 
approach  to  thermodynamics  where,  in  addition  to  the  usual  equations  for  conserva- 
tion of  mass  and  momentum,  an  equation  was  postulated  for  balance  of  entropy  and 
then  the  equation  for  balance  of  energy  was  required  to  be  satisfied  identically 
for  all  motions  and  all  temperature  distributions  after  the  elimination  of  the  ’ 
external  fields.  This  resulted  in  a procedure  for  obtaining  restrictions  on 


1. 


constitutive  equations  prior  to  any  consideration  of  a second  law  of  thermo- 
dynamics. The  particular  form  of  the  second  law  or  the  law  of  thermodynamic 
irreversibility  adopted  previously  (9a, b)  led  to  an  inequality  which  reflects 
the  fact  that  for  every  process  associated  with  a dissipative  material,  a part  , 

of  the  external  mechanical  work  done  is  always  nonnegative.  The  restriction  on 
the  heat  conduction  vector,  along  with  a restriction  on  the  internal  energy, 
were  discussed  previously  (£a)  on  the  basis  of  the  classical  concept  that  heat 
cannot  spontaneously  flow  from  a cold  to  a hot  part  of  the  body;  the  former 
restriction  was  obtained  from  consideration  of  the  heat  flux  response  in 
equilibrium  cases  for  which  heat  flow  is  steady,  while  the  latter  vas  deduced 
when  the  body  is  in  the  state  of  mechanical  equilibrium  with  the  temperature 
field  being  spatially  homogeneous . 

In  the  present  paper  we  extend  the  ideas  and  procedure  of  Green  and  Naghdi 
(£a,b)  to  mixtures  of  interacting  continua.  We  begin  in  $2  by  adopting  the  same 
conservation  laws  for  mass  and  momenta  as  those  employed  previously  by  Green  and  * 
Naghdi  (5)>  and  by  Craine,  Green  and  Naghdi  (2).  Further,  with  reference  to  thermal 
properties,  in  §3  we  admit  for  each  constituent  of  the  mixture  a different  field 
of  temperature  and  a different  field  of  entropy,  as  well  as  related  thermal 
fields;  and,  in  parallel  with  the  balances  of  momenta,  we  postulate  for  each 
constituent  a balance  of  entropy.  Next.,  we  recall  the  balance  of  energy  for  the 
mixture  aB  a whole  and,  after  elimination  of  external  body  forces  and  external 
entropy  supplies,  regard  the  resulting  equation  as  an  identity  to  be  satisfied 
for  all  thermo-mechanical  processes.  This  procedure  should  be  contrasted  with 

¥ 

existing  methods  for  mixtures  in  which  use  is  made  of  balances  of  energy  (rather 
than  a single  energy  equation  for  the  mixture  as  a whole)  anrt  in  general  leads 
to  different  constitutive  results.  For  example,  in  the  case  of  a mixture  of 
ideal  fluids  with  different  constituent  temperature  discussed  in  §4,  our  procedure 


2. 


yields  constitutive  relations  for  the  partial  pressures  and  diffusive  forces 
which  differ  from  those  found  by  Craine,  Green  and  Naghdi  (2),  except  when 
all  constituents  have  a common  temperature.  Inequalities  resulting  from  an 
appropriate  form  of  a second  law  for  mixtures,  as  well  as  restrictions  on  the 
heat  conduction  vectors  and  internal  energies,  are  discussed  in  §§3,5,6. 

We  return  to  a further  discussion  of  inviscid  fluids  when  the  mass  of  each 
constituent  is  separately  conserved  in  §6  and  also  examine  briefly  the  implica- 
tion of  the  second  law  and  restrictions  on  the  heat  conduction  vectors.  The 
next  section  (§7)  is  concerned  with  some  aspects  of  the  flow  of  an  incompressible 
Newtonian  viscous  fluid  through  a porous  solid  (taken  for  simplicity  as  rigid) 
whose  temperature  may  be  different  from  that  of  the  fluid.  Thermodynamic  restric- 
tions by  the  procedure  of  §3  are  obtained  for  various  coefficients  of  the  appropri- 
ate linearized  constitutive  equations,  representing  small  departures  from  an 
equilibrium  state  at  which  the  fluid  and  the  solid  have  the  same  constant 
temperature.  Also  included  in  the  development  of  §7  is  a brief  statement  of  a 
uniqueness  theorem  for  the  initial  boundary- value  problem  for  the  determination 
of  the  temperature  fields  in  the  fluid-rigid  solid;  and,  by  way  of  illustration, 
a discussion  of  heat  conduction  in  an  infinite  stationary  rigid  porous  solid 
containing  an  incompressible  viscous  fluid  at  rest.  We  conclude  the  present 
paper  with  some  additional  remarks  in  §8  concerning  certain  features  of  the 
present  and  the  previous  procedures  for  studying  the  thermodynamics  of  mixtures 
of  interacting  continue  and  also  comment  on  the  implU;  tion  of  an  entropy 
inequality  of  the  Claus ius-Duhem  type  for  .sr tores  as  cont.«. .«sted  with  the 
thermodynamic  restrictions  resulting  from  tht  procedure  proposed  here,  «s  well 
as  from  our  present  interpretation  of  the  second  . of  thermodynamics . 
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2. 


Mechanics  of  mixtures:  Summary  of  basic  equations. 


Consider  a mixture  consisting  of  v finite  constituents  C (a=  1,2,. . . ,v), 

Of 

each  regarded  as  a continuum  with  material  points  X°,  and  identify  the  material 
point  XQ  of  the  constituent  C with  its  position  Xa  in  some  fixed  reference 

Of  ~ 

conf iguration.  We  suppose  that  the  whole  of  some  region  of  Euclidean  space  ftQ 
is  simultaneously  occupied  by  all  constituents  in  their  reference  configurations. 
A motion  of  C is  defined  by  a sufficiently  smooth  vector  function  \a  which 

Of  ~ 

assigns  position  x® » xa(Xa,t)  to  each  material  point  Xa  at  each  instant  of 
time  t.  In  their  present  configurations  at  time  t,  all  constituents  Ca 
simultaneously  occupy  the  same  region  of  space  ft  which  varies  with  time  and 
which  is  bounded  by  a closed  surface  SR.  Similarly,  in  the  present  configura- 
tion, an  arbitrary  material  part  of  CQ  occupies  a portion  of  the  region  of 
space  ft,  which  we  denote  by  P (Cft)  bounded  by  a closed  surface  5P  . 

Of  “ Of 

Let  p^=  p^(X®,t)  be  the  mass  density  of  Cff  and  designate  its  velocity 
vector  at  time  t by 

where  Da(  )/Dt  denotes  material  time  derivative  holding  X®  fixed.  At  the 
place  x,  namely  at 

rs* 

x x (a  ■ 1,2,. ..,v)  , (2.2) 


Of 

v » 


Dt 


(2.1) 


we  define  the  density  p of  the  mixture,  the  mean  (or  the  barycentric)  velocity 
\ of  the  mixture  and  the  partial  diffusion  velocities,  respectively,  by 


v 

p - L 
or  1 


Pv 


v 

- E p v 
or -l  fl~ 


a 


(2.3) 


”inee  every  constituent  occupies  the  same  region  of  space  R at  each  instant  of 
time,  -it  tli1.;  pine,  x wo  have 


4. 


(2.4) 


v®  • n = v • n or  u®  • n = 0 (a=  1,2,. . . ,v)  on  aft  . 

where  n is  the  outward  unit  normal  to  3R. 

For  every  material  part  of  occupying  a region  p^  in  the  present  con- 
figuration, we  adopt  the  following  equations  for  conservation  of  mass,  momentum 
and  moment  of  momentum: 


p dv  •-  m dv  , 
J « of  J « « 


P 

a *p 

a 

, V®  dv  - f 

m 

v dv  - f ( 

a ~ Jp 

a 

a 

x®  x v®  dv  - I 
1 ~ ~ J 

m 

in  x®  x v dv 
a ~ ~ 

(2.5) 

(2.6) 


f [*“x(pb«-^)  + JLtf]dv+r  x®xt®da  . (2.7) 

J _ ri  ~ J ~ ~ 


In  (2.5)  to  (2.7),  m^  is  the  rate  at  which  mass  is  supplied  to  per  unit 
volume*  b®  is  the  external  body  force  per  unit  mass,  is  essentially  a 
measure  of  diffusive  force  per  unit  volume,  and  \®  is  the  internal  body  couple 
per  unit  volume. 

Making  the  usual  continuity  assumptions,  the  conservation  laws  (2.5)  to 
(2.',’’)  yield  the  local  forms 


if  + <iiv(PaI°f)  = ma  * 


div  T®+P  b“-u°  = m u®+p  f®  , t“*  T“n  , 

Of1**  Qt~  QJ'W  **  rw 


(2 

(2.9) 


This  may  include  a surface  contribution  across  dp  in  the  form  of  the  divergence 
of  a vector.  ® 

K 

The  detailed  development  of  the  conservation  laws  in  the  forms  (2.5)  to  (2.7) 
and  their  corresponding  local  forms  (2.8)  to  (2.10)  can  be  found  in  the  papers 
by  Green  and  Naghdi  ( 5 ) and  by  Craine,  Green  and  Naghdi  (2  ) . The  basic 

equations  (2.8)  to  (2.10)  are  equivalent  to,  but  not  the  same~as,  those  given 
by  Truesdell  and  Toupin  ( 10  ) . For  a discussion  of  equivalence  of  the  two 
systems  of  basic  equations,  see  Green  and  NBghdi  ( £ ) . Supplementary  discus- 
sion on  the  correlation  of  existing  forms  of  mixture’  theories  are  given  in  §6 
of  Green  and  Naghdi  ( 8 ) • 


lUVJi ' Vvtfv. ^ inV.vH A i v 1- • A.  '»- -i.. --  ..  , - 


(2.10) 


+ ip**  _ ip®  - o 

~ iV  /W  ^ 


where  all  variables  are  now  regarded  as  functions  of  x and  t,  div  stands  for 
the  divergence  operator  with  respect  to  the  place  x keeping  t fixed,  Ta  is  the 

__  r* 

T 

partial  stress  tensor  and  Ta  its  transpose,  n is  the  outward  unit  normal  to 

'■w  rw 

the  closed  surface  Sp^,  the  acceleration  vector  fa  is  given  by 


av*  av« 

f®  - + -r-  va 

~ at  ox  ~ 


(2.11) 


and  r°  is  the  skew  tensor  correponding  to  the  axial  vector  so  that 


r%  = \axu 


(2-12) 


for  every  vector  u. 

We  restrict  our  attention  here  to  a local  theory  of  mixtures  in  which  there 

is  no  mass  generation  within  the  mixture  and  in  which  the  internal  forces  and 

the  internal  couples  \a  do  not  give  rise  to  a resultant  force  and  a resultant 

couple  on  the  mixture  as  a whole.  Then,  for  material  volumes  of  each  constituent 

C which  coincide  at  time  t and  which  occupy  parts  P = p of  ft,  we  have 
o o 


L ”« dv  ■ 0 • A L 


o=l  P 


and  we  conclude  that 


a=l  *P 


dv  = 0 , £ \ dv  = 1 

o=l  Jp 
a 


£ m = 0 


Z a.  = 0 

_ r**  f+s 

o=l 


E \ = 0 
o=l  ~ ~ 


, e r = o 

0=1  ~ ~ 


(2.13) 
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3*  Thermal  properties,  Themomechanics  of  mixtures. 

We  consider  here  the  thermal  properties  of  the  mixture  and,  for  each 

constituent  C , introduce  first  the  absolute  temperature  at  each  material 

point  by  a scalar  field  0q  = 0^(X*,t)>O.  Along  with  the  temperature,  we  admit 

the  existence  of  an  external  rate  of  supply  of  heat  ra=r^(xa,t)  per  unit  mass 

and  an  external  rate  of  surface  supply  of  heat  - per  unit  area  acting  across 

dR.  Also,  we  assume  the  existence  of  an  internal  surface  flux  of  heat 

-h  *-h  (Xa,t  ; n)  per  unit  area  across  each  surface  dp  ; the  field  h , called 
a a~  ~ cr  a 

the  partial  heat  flux  and  measured  per  unit  area  per  unit  time,  assumes  the 

value  h on  dR.  We  now  define  the  ratio  of  heat  supply  r to  the  temperature  0 
a Of  Of 

as  s *e  (X®,t)  and  call  this  the  partial  external  rate  of  supply  of  entropy  per 
unit  mass.  Similarly,  we  define  the  ratios  of  and  to  the  temperature  0 , 
respectively,  as  the  partial  external  rate  of  Burfa.ce  supply  of  entropy  k^  per 
unit  area  of  dR  and  the  partial  internal  surface  flux  of  entropy  k =k  (Xa,t  ; n) 
per  unit  area  of  dp^.  These  definitions  may  conveniently  be  summarized  by 


0S  , h = 0 k ,h-0k 
q?  of  a a a a a a 


(3.1) 


In  addition,  for  each  constituent  C we  assume  the  existence  of  a scalar  field 

Of 

7]  = T|a(Xa,t)  per  unit  mass,  called  the  partial  specific  entropy , and  a 
partial  internal  rate  of  production  of  entropy  5^  = 5^(Xa,t)  per  unit  mass. 

The  contribution  of  the  latter  to  the  internal  rate  of  production  of  heat  is 
simply  QJ~>a  per  unit  mass. 

We  now  postulate  a balance  of  entropy  for  every  material  volume  of  Cq 
occupying  a part  Pq  in  the  present  configuration  and  write 

— f p T|  dv  - f p (s  + P ulv  - f k da  . (j 

dt  J , 'a  V a *a)  a 

a 1 of  a 

By  usual  procedures  it  can  be  shown  from  (3*2)  that  k is  linear  in  n,  i.e., 


7. 


k = * n j 

a ~ ~ 


(n  'j  ' 


where  p 


a is  called  the  partial  entropy  flux  vector.  Then,  from  (3-lL  and  (3-3) 


a 


, = 0 pa  . n and  we  may  define  the  partial  heat  flux  vector  q by 
a cfi'  ~ 


a 


a 


a “ e<£  • 

Under  suitable  continuity  assumptions  and  with  the  use  of  (3*3);  the  local, 
equation  at  the  place  x resulting  from  (3-2)  is 


(3.*0 


DaT] 


a 


j + m Tt  r:  p . , 

a Dt  of  a a a ot 


(s  + F ) - div  pa 


(3 .5) 


whefe 




Dt 


an 


dt 


'Of  , Of 
• + V 


(3.6) 


At  this  point  we  introduce  the  first  law  of  thermodynamics  or  the  principle 
of  balance  of  energy,  which  states  that  heat  and  mechanical  energy  are  equiva- 
lent and  together  they  are  conserved.  We  state  the  law  of  conservation  of 
energy  only  for  the  mixture  as  a whole.  However,  the  conservation  of  energy 
could  also  be  employed  for  each  constituent  provided  full  allowance  is  made 
for  interaction  energy  from  all  other  constituents,  but  we  do  not  consider  this 


heret  Thus,  with  reference  to  the  present  configuration  and  for  material  volumes 
which  coincide  at  time  t and  occupy  parts  P^  = P of  ft,  the  balance  ot  energy  for 
the  mixture  may  be  stated  in  the  form 


This  is  in  contrast  with  previous  work  on  many  temperatures  by  Green  and 
Naghdi  (4)  and  by  Craine,  Green  and  Naghdi  (2),  where  it  was  necessary  to  have 
energy  balances  for  each  constituent  in  order  to  have  sufficient  equations 
for  the  temperature  fields. 
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(3.7) 


E 7T  \ (2  P v°  . v®  + p < )dv 

a=i  dt  Jp  ~ at  at 

a v 

= E f ( p b®  . v®  + p r )dv  + £ f (t®  • v®  - h )da  , 

, J — of-  ~ at  a , ~ ~ a 

or=l  P a=l  "dp 

a a 

where  *aa *^(X®,t)  is  the  partial  specific  internal  energy  of  the  constituent 
Ca>  With  the  help  of  (2.8),  (2.9) > (3.1) » (3«*0  and  (3-0)  end  under  appropriate 
continuity  assumptions,  the  local  energy  equation  at  the  place  x resulting  from 

/w 

(3.7)  is 


v D®«  D®TL 

0<.‘tT)+"a<,«-9A)1 


+ E Cp.01  • v®  + T®  . grad  v®  + f m u®  • u®] 

Qf-1 


- 2 tpaSA+E  4 grad  6J#0  » (3.8) 

Qf=l  B 01  “ W 

where  all  variables  in  (3*8)  are  regarded  as  functions  of  x end  t and  where 
grad  stands  for  the  gradient  operator  with  respect  to  x keeping  t fixed. 
Introducing  the  specific  Helmholtz  free  energy  “ ^(X^t)  by 


VV'A  ■ 


(3.9) 


the  energy  equation  (3*8)  may  be  written  in  the  alternative  form 
v D®i|f  D®Q 


+ E [u®  • v®+  T®  • grad  va  + £ m u®  . u®l 

K#  r*  w rv  e ryw  rw 

a=l 

v 

“ 2 9 +P  • grad  9 ] = 0 . 

1 at  at  at  ^ at 


(3.10) 


a*l 

For  a given  mixture  of  v constituents  having  reference  mass  densities 

p (X®),  she  J.o”il  field  equations  resulting  from  the  integral  forms  of  the 
ao  '■w 

conservation  laws  involve  a set  of  12 v functions.  These  consist  of  the 

d- 'ormation  functions  X®  and  the  temperatures  9 , i.e., 

~ 0 


9. 


nam'd y 


and  the  various  mechanical  and  thermal  fields, 


tViV.itVV  VE“5 


,3-12) 


and 


{ba,s  } ( 3 • 1 3 J 

L~  a 

The  mass  densities  p are  not  included  in  (3.1-2)  since,  given  (3-11),  they  can  be- 
ar 

calculated  from  (2.8).  We  assume  that  the  fields  (3.12)  are  specified  by 
constitutive  equations  which  may  depend  on  the  variables  (3.11),  their  space 
and  time  derivatives,  as  well  as  the  whole  history  of  deformation  and  temperatures. 
We  then  adopt  the  following  procedure  in  utilizing  the  conservation  laws: 

(1)  The  field  equations  are  assumed  to  hold  for  an  arbitrary  choice  of 
the  functions  (3. 11)  including,  if  required,  an  arbitrary  choice  of  the  space 
and  time  derivatives  of  these  functions; 

(2)  The  fields  (3.12)  are  calculated,  from  their  respective  constitutive 
equations ; 

(3)  The  values  of  ba  and  s can  then  be  found  from  the  balances  of 

~ or 

momenta  (2.9)  and  the  balances  of  entropy  (3«5)> 

(4)  The  conditions  (2.13) > the  equations  (2.10)  which  are  deduced  from 
the  balances  of  moment  of  momenta  and  the  energy  equation  {3. 10}  will  be 
regarded  as  identities  for  every  choice  of  (3. 11).  These  wilD  place  restric- 
tions on  the  constitutive  equations. 

In  the  current  literature  on  continuum  mechanics  sometimes  certain  terminology 
such  as  a thermodynamic  process  (or  simply  a process ) and  an  admissible  process 
are  employed.  In  the  context  of  the  present  paper,  the  former  refers  to  the 
specification  of  12v  functions  (3.11)  to  (3. 17,)  such  that  they  are  compatible 


with  the  field  equations  resulting  from  the  conservation  laws  (2.5),  (2.6), 

(2.7)  and  (3.2).  In  order  to  specify  a thermodynamic  process  it  will  suffice 
to  prescribe  the  lOv  functions  (3«ll)  and  (3*12)  and  the  remaining  2v  functions 
(3.13)  are  determined  by  the  field  equations  resulting  from  (2.6)  and  (3*2), 
which  corresponds  to  (3)  above.  A thermodynamic  process  is  said  to  be  admissible 
in  the  mixture  if  it  is  compatible  with  the  constitutive  assumption  at  each  point 
of  the  mixture  for  all  times,  and  this  corresponds  to  (1)  and  (2)  in  the  above 
procedure.  The  last  step  (4)  then  places  restrictions  on  admissible  processes. 

Some  additional  remarks  are  necessary  concerning  the  entropy  equations 
(’.5).  Suppose  that  the  external  fields  (3*13)  have  prescribed  values  and 
that  the  constitutive  equations  have  been  specified  for  the  set  (3.12)  in 
terms  of  (3.11)  and  their  histories  by 


m = 
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T«,£r 


a *0 
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a 

P * P 


, c “ ♦ + 8 1) 

' a fa  a 1 a 


(3-14) 


Then,  the  requirement  that  the  energy  equation  (3. 10)  be  identically  satisfied 
leads  to  the  following  identity 


ft  D*8- 


♦ * v°  + T®  • grud  ,«i  ’ 

~ r%J  w ^ Q v**  ^ 


o-l 

v A (V 


■ E tp«V«+£  * grad  e«]  = 0 
0=1 


(3.15) 


relating  the  response  functions  in  (3 -14).  Also,  the  deformation  functions 
and  the  temperatures  must  satisfy  the  differentia]  equations  which  result  from 
substitution  of  (3.14)  into  (2.g)1  and  (3.5).  Now,  let  the  constitutive 
equations  be  changed  to 


U. 


m = m , 
or  a 


T1  - H + f 
'a  a a 


’Of 
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a Aa  ,o f {Vx 

T = T , 

U = P > x.  = X 

/V  cw  /w 
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Of  Of 
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- 0 f , 

* rv  rj 

g - ik  + 0 T)  , 

(3.16) 


where  f are  arbitrary  scalar  functions  of  the  variables  (3.1l)>  their  space 
and  time  derivatives  and  functionals  of  their  past  histories.  The  partial 
stress  tensors,  the  diffusive  forces  and  couples,  the  entropy  flux  vectors 
(and  hence  the  heat  flux  vectors)  and  the  partial  specific  internal  energies 
(3.14)^  are  unchanged  and  the  energy  equation  (3*15)  which  imposes  restrictions 
on  the  response  functionals  in  (3.14)  is  also  unchanged.  If  now  (3.16)^  to 
(3»l6)g  are  substituted  into  (2.9)^  and  (3.5),  we  obtain  the  same  differential 
equations  for  the  variables  (3«H)  as  those  resulting  from  the  set  (3*14). 

Also  the  boundary  conditions  imposed  on  Xa  or  T°  and  on  0®  or  p®,  together 
with  appropriate  initial  conditions  on  X®  and  0 , are  unaltered.  Once  the 

~ Of 

solutions  for  the  deformations  and  temperatures  have  been  determined,  there 

remains  some  arbitrariness  in  D ,5  ,4  to  the  extent  of  the  functionals  f in 

■of  of  of  a 

(3.16),  even  though  the  differential  equations  and  the  boundary  and  the  initial 
conditions  do  not  involve  f^.  An  examination  of  the  balance  of  entropy  for  the 
ath  constituent  as  given  by  (3-5)  reveals  that  f represents  an  internal  genera- 
tion of  entropy  which  is  accounted  in  (3.5)  by  an  internal  rate  of  production  of 
entropy  5 . Moreover,  this  internal  rate  of  production  of  entropy  can  be  sustained 

in  any  part  P without  change  in  the  partial  specific  internal  energy  c and  with 
a a 

no  supply  of  entropy  from  external  sources.  This  is  a form  of  internal  generation 

of  entropy  which  seems  to  be  physically  unrealistic.  In  order  to  exclude  this, 

we  remove  the  arbitrariness  in  f by  setting 

a 

D°f 

*«  -Dt2  * V«  ■ 0 ’ "A b.  * ’ 
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(3.17) 


where  F is  the  partial  deformation  gradient  defined  by  F *dxa/3X0  and  f 

I i\Gt  rv  ' i N/  I 


~o 


oo 


is  an  arbitrary  function  of  Xa.  Then,  is  determined  uniquely  and  7)  is  only 


arbitrary  to  the  extent  of  an  additive  function  f given  by  (3.17).  The  function 

o ' 

fao(X°)>  and  hence  f^*  is  determined  by  specifying  a value  for  T|  in  some 
reference  state. 


The  cace  in  which  all  constituents  are  at  a common  temperature  e is  of 
particular  interest  and  may  be  regarded  as  a special  case  of  the  above  theory 
in  which  all  constituent  temperatures  are  constrained  to  be  equal,  i.e., 


ea  “ ? (a  = 1,2,. ..,v) 


(3.18) 


A constraint  of  this  kind  is  possible  only  if  v-1  of  the  external  rates  of 
supply  of  entropy  are  chosen  so  as  to  balance  v-1  of  the  entropy  equations  (3*5) 
while  the  total  external  rate  of  supply  of  entropy 


ps  * E p s 


a=l 


a a 


(3.19) 


remains  at  our  disposal.  Equations  (3.2)  are  then  replaced  by  a single  equation 
for  the  mixture  as  a whole  in  the  form 


i « L dv  • L p(3j?)dv  ■ J9P  k 


da 


(3-20) 
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where 


p?  = E p 5 , k-Ek-p-n.p-JEp 

o=l  a a o=l  a ~ ~ ~ o=l  ‘ 
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(3.21) 


and  where  the  total  heat  flux  and  the  tote  L heat,  flux  vector  are  given  by 

v 


h ■■■■  & ■ a > i- 1 -- 
0=1 


(3.22) 


The  local  field  equation  which  follows  from  (3*20)  r<’uds  as 


13. 
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(3.23) 
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and  the  energy  equation  (3. 10)  now  becomes 
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(3. 2k) 


Before  considering  what  restrictions  arise  from  some  form  of  second  law 

of  thermodynamics,  we  discuss  in  the  next  section  a special  case  of  the  above 

thermo -mechanical  theory  for  a mixture  of  inviscid  fluids. 

For  later  use,  we  need  to  record  the  expressions  for  the  external  mechanical 

work  and  the  external  heat  supplied  to  the  whole  mixture  occupying  the  region  of 

space  ft  during  the  time  interval  t^StStg.  First,  however,  we  observe  that  in 

the  case  of  the  mixture  of  inviscid  fluids  discussed  in  §4  the  response  functions 

for  depend  only  on  the  constituent  densities  p^  and  the  constituent 

tenjperatures  9 but  are  independent  of  their  rates,  the  velocities  v^  and  the 
P 

« 

temperature  gradients  grad  9 . An  elastic  mixture  of  this  kind  will  be  regarded 

P 

as  nondissipative  in  a sense  that  will  be  made  precise  later;  and,  in  conjunction 
with  an  expression  for  the  external  mechanical  work  supplied  to  ft,  will  be  used 
as  a basis  for  establishing  in  §5  an  inequality  representing  the  second  law  of 
thermodynamics  for  dissipative  materials.  Keeping  this  background  in  mind,  we 
assume  that  the  constitutive  response  functions  for  *a>'Ha  include  also  dependence 
on  the  set  of  variables  v“,9  ,grad  @ and  their  higher  space  and  time  derivatives 

~ Q Of 

and  refer  to  this  set  collectively  as  \j.  Further,  let  denote  the  respective 

values  of  when  the  set  V is  put  equal  to  zero  in  the  response  functions. 

Thus,  for  example, 


These  definitions  of  © r , T\ ' do  not  exclude  their  dependence  on  the  past 


histories . 


Of  'Of 
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V = (v® ,D®0 ^Dt,grad  0^, . . . ) , 


where  the  dots  in  (3.25)2  refer  to  the  higher  space  and  time  derivatives  of 


va>Da0o/Dt,grad  9^.  Then,  with  the  help  of  (2.9)  and  the  integral  of  (3.8)  with 
respect  to  time,  we  obtain 
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and 
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where 
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4.  Mixture  of  lnvlscid  fluids. 


Consider  a mixture  of  inviscid  fluids  with  v constituents.  We  define  the 
properties  of  each  constituent  by  the  constitutive  assumptions  that  the  set  of 
variables  (3*12)  are  functions  of 


Pp»  8p,  grad  pp,  grad  (p  = 1,2, . . . ,v)  . 


(4.1) 


As  indicated  in  $3  [following  ( 3 • 13 ) ] # the  equations  (2.10),  the  conditions 
(2.13)  and  the  energy  equations  (3*10)  are  regarded  here  as  identities  to  be 
satisfied  by  all  motions  and  all  temperature  distributions . Using  this 
procedure  and  making  use  of  (2.8),  the  application  of  standard  techniques  leads 
to 
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(4.5) 


(4.6) 


where  I is  the  identity  tensor  and  in  writing  (4.6)  we  have  also  recalled  (3.28)^. 


In  view  of  (4.6),  we  put 


w. = - ia  • ^ 9a + w;  . i P0eos;  - o 


(4.7) 


Then,  from  (3«5)  and  (4.7)  we  have 
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where  the  rate  of  supply  of  heat  for  the  mixture  and  the  total  heat  flux  vector 
are  defined  by 

v v 


Pr  “ E P * » % “ 2 a • 

o-l  a * «-l 

For  a mixture  of  inviscid  fluids  characterized  by  (4.4)  to  (4.8),  with  the 
help  of  (4.8)  and  the  results  of  §§2,3>  we  record  below  the  egressions  for  the 
mixture  as  a whole  representing  (i)  external  work  by  body  and  surface  forces  on 
the  mixture  and  (ii)  external  supply  of  energy  to  the  mixture  arising  from  the 
rate  of  supply  of  heat  and  the  total  surface  flux  of  heat,  both  over  a finite 
time  Interval  t-^Stetgj 


(4.9) 
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and 
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In  the  above  discussion  of  a mixture  of  inviscid  fluids  the  response  func- 
tions are  assumed  to  depend  only  on  the  set  of  variabl  cs  (4.1)  which  excludes 
the  velocity  vectors  va.  An  extension  of  the  discussion  to  include  dependence 

rw 

on  the  velocity  vectors  is  carr  .ed  out  in  §6.  In  the  meantime,  we  observe  that 
the  development  of  the  present  section  can  be  extended  to  include  a mixture  of 
elastic  solids,  or  elastic  solids  and  inviscid  fluids.  In  general,  the 


17. 


constitutive  assumption  for  a mixture  of  elastic  materials  requires  that  the 
variables  (3*12)  be  functions  of  the  temperature,  temperature  gradients, 
deformation  gradients  and  second  deformation  gradients.  With  the  help  of 

(2.10) ,  (2.13)  and  (3.10),  it  follows  that  the  response  functions  for 

depend  only  on  the  temperatures  and  deformation  gradients.  With  these  resultB 
and  recalling  the  notations  of  §3,  we  may  put  \ji^  = f sind  it  can  then 

be  shown  that  the  expressions  for  the  external  mechanical  work  to  and  the  external 
heat  H supplied  to  the  mixture  again  reduce  to  the  forms  (4.10)  and  (4.11), 
respectively.  We  make  use  of  these  results  in  §5. 

In  the  case  in  which  all  constituent  temperatures  are  equal,  instead  of 

(3.10) ,  the  energy  equation  is  given  by  (3*24)  and  the  results  (4.2)  to  (4.6) 
are  replaced  by 
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Also,  instead  of  (4.8)  , we  now  have 
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18. 


For  a mixture  with  many  constituent  temperatures,  Craine,  Qreen  and 
Naghdi  ( 2 ) have  discussed  more  general  constitutive  equations  than  is 
implied  by  the  dependence  of  (3.12)  on  the  variables  (4.1).  The  detailed 
development  of  their  work  on  a mixture  of  viscous  elastic  materials  is  such 
that  direct  comparison  of  their  results  with  (4.2)  to  (4.5)  is  difficult. 
However,  if  the  energy  equation  (3.10)  is  not  regarded  as  an  identity  and  the 
entropy  inequality  employed  by  Craine  et  al.  ( 2 ) is  used  to  impose  restric- 
tions on  the  forms  of  Ta,i)t  ,1)  ,ua,X.a  as  functions  of  the  variables  (4.1),  then 
the  following  results  emerge: 
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The  above  resullr  differ  from  those  by  (4.2)  to  (4.5),  although  the 

differences  disappear  when  *1j.  the  temperatures  ure  equal.  t+  r ho  ala  tn-. 
emphasized  that  the  formulae  (’i.2)  to  (4.5)  have  been  obtained  without  any 
appeal  to  a second  law  of  thermodynamics,  in  contrast  to  the  results  (4.18) 
to  (4.22),  which  are  obtained  here  in  the  spirit  of  euili.r  procedures  with 
the  use  of  an  entropy  inequality. 


(4.18) 

(4.19) 

(4.20) 


(4.21) 
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19. 


i 
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vectors  and  internal  energies. 

Previously  [see  Green  and  Naghdi  (9)>  §4],  we  examined  the  nature  of 

thermodynamic  irreversibility  arising  from  a mathematical  interpretation  of  a 

statement  of  the  second  law  of  thermodynamics , namely  that  "it  is  impossible 

completely  to  reverse  a process  in  which  energy  is  transformed  into  heat  by 

friction."  nils  led  us  to  postulate  an  inequality  which  is  different  from 

the  Claus ius-Duhem  inequality;  some  of  the  undesirable  features  of  the  latter 

inequality  are  discussed  in  $7  of  Green  end  Naghdi  (£). 

Here  we  follow  a procedure  for  mixtures  which  is  similar  to  that  used 

previously  for  single  phase  continue.  First,  we  observe  that  the  expressions 

(3.26 )1,  when  evaluated  for  a given  process,  may  be  either  positive  or  negative 

depending  on  whether  the  external  work  is  supplied  to,  or  is  withdrawn  from,  R. 

This  external  work  is  also  represented  by  (3.26 )2  in  terms  of  both  thermal  and 

mechanical  variables  but  not  every  term  in  (3.26)g  need  necessarily  be  positive 

(zero  or  negative),  even  though  the  external  work  may  be  positive  (zero  or 

negative).  Thus,  following  a discussion  similar  to  that  in  (2,10) , with 

reference  to  a dissipative  mixture  we  assume  that  in  every  admissible  process 

a part  of  the  external  work  done  is  always  nonnegatlve.  Then,  in  a process  end 

its  reverse  process  at  the  end  of  which  the  mixture  has  returned  to  its  same 
ft 

state,  some  of  the  work  done  is  always  transformed  into  heat.  Hence,  we  write 

to  = u^  + ih,  , . 

The  above  also  implies  that 

Id  £ I&l 

1. .v.'  it  mounded  below,  the  bound  depending  on  the  process. 

f’O. 


In  order  to  identify  the  two  parts  \a,  and  we  note  from  (4.6)  and  (3.26) 

that  in  the  case  of  a mixture  of  inviacid  fluids  or  elastic  materials  discussed 

in  §4  the  expression  for  the  external  mechanical  work  supplied  to  R reduces  to 

(4.10)!  We  now  regard  an  elastic  mixture  of  this  kind  as  nondissipative  in  the 

sense  that  no  restriction  is  placed  on  the  external  mechanical  work  supplied  to 

v 

R and  identify  lo,  with  the  right-hand  side  of  (3.26),,  after  setting  £ p w =0. 

” 2 a-1  “ “ 

Keeping  this  in  mind*  we  rewrite  the  last  inequality  as 

v c - . I*a  ;t0i.  v a(pji') 


(5.1) 


lb  * £ f p (£  va  • va  + « )dv  2-f  2f  £ 8 [ — £ + div(  p T\'va)  )dv  dt 

and  assume  that  (5.1)  holds  for  every  thermo-mechanical  process.  The  combination 


of  (3«26)2  and  the  assumption  (5*1)  yields 


nZ  p w dv  dt  £ 0 

a 01  a 


't^R  Od 


(5-2) 


or*l 


for  all  times  t^,tg  (t^aitg).  Since  t^tg  are  arbitrary  and  £ p^w 
already  been  assumed  to  be  continuous,  it  follows  that 

r p w dv  £ o 


has 


v 
E 

«-l  "R 


(5.3) 


for  all  thermo-mechanical  processes.  With  the  help  of  (3.27)  and  (5*3) , we 
also  have 


v 

was 


ri. 


0=1  Jt1JR  0=1  * 


v a(p  tT) 

2 9 J — + div ( p J\V* ) }dv  dt 


3t  ■ ro  'ar~ 

so  that  the  external  heat  supplied  to  the  whole  mixture  is  bounded  above  in 
every  process.  Alternatively,  some  of  the  heat  supplied  to  the  mixture  in 
every  process  is  always  nonpositive. 


(5.4) 


^It  may  be  recalled  that  the  expressions  (4.10)  and  (4.11)  hold  also  for  a 
mixture  of  elastic  solids  (or  fluid-solid)  in  which  dependence  on  the  velocity 
vectors  is  excluded,  as  was  noted  in  the  paragraph  following  (4.11). 


21. 


The  development  of  this  section  began  with  an  examination  of  the  expression 
for  the  supply  of  external  mechanical  work  to  the  whole  mixture  occupying  the 
region  R at  time  t and  this  led  us  to  postulate  the  inequality  (5.1)  from  which 
the  thermodynamic  condition  (5*3)  was  deduced.  It  is  also  possible  to  obtain  a 
similar  condition  for  every  part  of  the  mixture.  For  the  latter,  instead  of 
considering  the  total  mechanical  work  during  a given  time  Interval,  we  consider 
the  rate  at  which  mechanical  work  is  supplied  to  a part  of  the  mixture.  Thus, 
with  reference  to  the  present  configuration  and  for  material  volumes  which 
coincide  at  time  t and  occupy  parts  p^=  p of  R,  the  rate  W(t)  at  which  external 
work  is  supplied  to  a part  p of  the  mixture  is 

W(t)  <*  E f p b®  • vadv  + £ f ta*v“da  . (5*5) 

awl  up  * a*l  "dp 

a 

Using  results  similar  to  (3.26)  and  (3.28)  but  with  R replaced  with  P and  the 
limit  of  integration  t0  replaced  with  t,  the  last  expression  can  be  put  in  the  form 


v 

W(t)  = z 
a=>l 


£ J <* 


P v 


a 


va+  p c 
~ Ka  a 


)dv 


v p 

+ z J C-e  C 

awl  J p * 


St 


+ div(p  h'vQ')]  + p w ]dv 
a a a 


(5.6) 


By  an  argument  similar  to  that  which  led  to  (5.1),  for  each  part  P of  the 
mixture  we  write 


W(t)  5 £ 
awl 


•rr  I (|  p va  . v°  + p e )dv 

it  ‘ ra~  ~ Ka  a 


v r 

■ E,  j„  9«[— + “''(p.V. >]iv 

o/=l  r 


(5-7) 


and  assume  that  (y.,)  i ;■.> I i.-.  for  every  thermo-mechunical  process.  The  combina- 
tion of  (5*6)  and  the  assumption  (5.Y)  yields 


v 

E 

a=l 


p w dv  £ 0 
a a 


(5-8) 


22. 


for  every  part  p.  From  this  last  inequality  follows  the  local  form 

v 

£ P w 2 0 , (5.9) 

, a 1 a 
Ofsl 

which  must  hold  for  all  thermo-mechanical  processes. 

We  supplement  the  above  interpretation  of  the  second  law  with  two  further 
inequalities  which  stem  from  the  statement  that  "heat  cannot  pass  spontaneously 
from  a body  of  lower  temperature  to  a body  of  higher  temperature."  To  this  end, 
we  consider  first  the  heat  flux  response  in  equilibrium  cases  for  which  the  heat 
flow  is  steady.  By  equilibrium  we  mean  that 


s“-«  • • 8“'9  - V f°r  a11 1 * <5-10> 

where  F®  is  the  deformation  gradient  for  the  constituent  C » the  tensor  F®, 

~ 9 ~0 

the  common  temperature  9 and  the  densities  as  well  aB  all  other  relevant 

functions,  are  independent  of  t but  may  depend  on  x.  Then,  m = 0 and  the 

~ a 

conservation  equations  (3. 5)  and  (3*8)  become 

v 

Po/B«+Sj  = dlV  P°  > S (Pj«e  + p0(-  grad  6)  = 0 ‘ (5-11) 

ur  u Or  **  ® ^ 

Also,  since 


pr  * E or  = £ p0e  = 0 £ p s 

" in  fit  ry  r ™ 


9=1  ® ® 9,1 


% - Z <f  = 2 9 P' 

9=1  9=1 

it  follows  from  (5.H)  that 


999 
9 


9=1 
v 

0 £ p 

1 ~ 

9=1 


9 9 


9 


(5.12) 


pr  = div  q . (5.13) 

The  equation  (5<13)>  which  is  obtained  hi.:re  for  a mixture  in  equilibrium, 
has  the  same  form  as  the  corresponding  result  for  a single  phase  continuum; 
and  thus,  for  equilibrium  cases  under  discussion,  wc  again  adopt  the  classical 
heat  conduction  inequality 
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- q • grad  0 £ 0 


(5.14) 


for  all  time- independent  temperature  fields  (tj.lO)^.  We  note  that  the  inequality 
(5.14)  can  be  written  in  the  alternative  form 


v 

- £ P°  • grad  0 £ 0 , (5-15) 

«=1  ~ 


in  view  of  (5.12)  and  since  0>O. 

In  the  remainder  of  this  section,  we  suppose  that  the  mixture  is  at  rest 

with  va=0  for  all  time  and  with  the  densities  p and  with  the  deformation 
~ ~ ct 

gradients  dX°VdXa  everywhere  constant  for  all  time.  Then  m =0  and  the  total 

rv  & 

V 

mass  density  p = t p , which  is  positive,  is  constant.  In  addition,  we 
o=l  a 

restrict  the  temperature  fields  to  be  spatially  homogeneous  and  equal  so  that 

ft  =0(t).  The  specific  internal  energies  e which  will  depend  on  p , }sxa/'$Xa, 
a 0 a ~ 

6 and  its  history  will  be  functions  of  t, independent  of  x.  Moreover, 
recalling  (4.9),  from  combination  of  (3*5)  and  (3.8)  we  have 

a v 

pr-div  q = ■&;  (pe)  , pe  = pfflea  , (5.16) 

where  e is  a function  of  t independent  of  x.  Since  v0  = 0 everywhere,  no 
mechanical  work  is  supplied  to  the  whole  mixture  occupying  the  region  R at 
time  t.  Then,  by  integration  of  (5.16^  with  respect  to  t,  the  heat  supplied 
to  R during  the  time  interval  t^  * t 5 is 


P 


(5.17) 


We  now  suppose  that  the  mixture  has  been  in  thermal  equilibrium  during  some 
interval  of  time  up  to  the  time  t^  with  constant  temperature  9-^  and  constant 
Internal  energy  7.  We  then  assume  that  whenever  heat  is  supplied  to  the  mixture 
according  to  (5. 17),  the  temperature  0(t)  throughout  R will  be  increased.  Hence, 
we  impose  the  condition  that 


2k. 


2 > 0 whenever  W > 0 . (5*18) 

i 

rws  from  (5*17)  and  (5*18)  that 
> 0 whenever  c(t)  -«  > 0 (5«19) 


6.  Further  discussion  of  inviscid  fluids. 


In  order  to  illustrate  the  nature  of  the  preceding  thermodynamical 
developments,  we  return  to  a mixture  of  inviscid  fluids  which  is  partially 
discussed  in  §4.  For  simplicity,  however,  we  limit  attention  to  a mixture  of 
two  fluid  constituents  whose  masses  are  separately  conserved.  Thus,  we  write 


DXp 


sr + pi div  i1  * 0 


D% 

”2  ■ “ + H 


div  v o 0 


(6.1) 


and,  in  view  of  (2.13)  , we  adopt  the  notations 

<-  >3 


(6.2) 


Constitutive  equations  for  such  fluids  usually  include  the  dependence  on  the 
velocity  vectors  but  these  were  omitted  in  the  discussion  of  section  4.  We 
now  assume  that  the  set  of  variables 

£ >I2 * * * *>  V 

are  functions  of 


px»  P2>  erad  P!>  STS'1  P2»  91»  ®2»  0rad  Krad  &2  (6.4) 

1 2 

and  v ,v  . In  view  of  invariance  conditions  under  superposed  rigid  body  motions, 

r*  rs* 

1 2 

however,  the  variables  v ,v  must  be  replaced  by 


1 

v 


= a 


(6.5) 


In  accordance  with  the  procedure  of  §3,  since  (2.10),  (2.13)  and  the  energy 

equation  (3.10)  must  be  satisfied  as  identities  for  all  motions  and  all  tempera- 

ture distributions,  we  then  obtain 
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*1  “ * *2  " *2^pl»%»ei»®2*!^ 

P*  " Pl*l+p2*2  » * “ ^P1>P2>e1>92)  > 


»1%  ■ ■ » 

A 

*♦1 


S*, 


A 

3*1 


2?--<*%sH&-Pi£a  • 


fc  * £+ Pi  gr*4  •fe'Ps5^«r,d  ”i 

3*1  3*P 

+ pi  ae2  s™*  ®2  ’ p2  ae1  8rRd  0i  * 


1 2 — 

P1?191+  p2?2^  +E  * 9l  +E  * grad  0p  - J*  • a - 0 , 


(6.6) 


(6.7) 


(6.8) 


(6.9) 


where  j*  Is  a function  of  the  variables  ( 6.4 ) and  (6.5).  Apert  from  invariance 

considerations  no  further  reduction  in  the  form  of  these  results  is  possible 

without  the  use  of  the  thermodynamic  restrictions  (5.3)  or  (5.9)  and  (5.15 ) or 

without  employing  more  restrictive  constitutive  assumptions.  Here  we  uec  the 

inequality  (5*9)  with  w given  by  (3.28).  From  this  inequality,  we  may  deduce 

o 

the  result 


1 t 0 (6.10) 

for  all  thermo -mechanical  processes.  Now  put 

fc"fc/  + ? » 

where  y/  is  the  part  of  the  response  function  for  ^ when  a is  set  equal  to 
zero.  Then,  assuming  that  £ is  a continuous  function  of  its  arguments,  from 
(6.10)  we  deduce  the  results 
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(6.11) 


n'  ■ 2,  » • a. * 0 • 

Al  ho 

Also  some  restrictions  on  the  forms  of  the  response  functions  p ,j?  may  be 
found  with  the  help  of  (5.15) > but  we  leave  the  discussion  at  this  point. 


7.  Flow  of  a fluid  through  a rigid  solid. 

In  many  problems  involving  the  flow  of  a fluid  through  a porous  solid,  the 
latter  is  often  regarded  as  a rigid,  rather  than  deformable,  solid;  and,  in 
addition,  both  media  are  assumed  to  be  at  the  same  temperature.  Our  main 
interest  here  is  to  explore  a further  application  of  the  theory  of  earlier 
sections  with  reference  to  the  flow  of  a homogeneous  incompressible  viscous 
fluid  through  a porous  solid  when  the  fluid  and  the  solid  are  at  different 
constituent  temperatures;  and,  in  order  to  avoid  undue  complications,  we  con- 
sider the  solid  to  be  stationary,  homogeneous  and  rigid. 

Let  the  incompressible  viscous  fluid,  with  its  mass  element  conserved, 
be  referred  to  as  the  constituent  and  let  the  stationary,  homogeneous  and 
rigid  solid  be  identified  as  the  constituent  v2*  Then,  ^ and  pg  are  constants, 


mi  = 0 , mg 


1 2 

m„  * 0 , div  v - 0 , v -0 


(7-1) 


and  by  (2.13 )g  3 4 we  also  have 


x2--*1  , r2  « - r1 


(7.2) 


The  energy  equation  (3. 10)  in  this  case  becomes 


D1*  D191  dtP  d9? 

• f>i<  ^ + % Tr>  • '2(ir + \ 


■*  & - l er»d  l -fAh  - i>2?202 


1 2 

- 53  • grad  01  - £ • grad  0^  0 


(7.3) 


Because  of  the  incompressibility  condition  (7-l)3»  the  stress  tensor  T is 
undetermined  to  the  extent  of  an  additive  hydrostutic  pressure  p so  that 


TL  = - pi  + T , 


(7.4) 
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where  p is  an  arbitrary  scalar  function  of  x,t  and  T is  to  be  specified  by  a 
constitutive  equation. 

For  the  two  constituents  under  consideration,  we  introduce  the  constitutive 

assumptions  that 

rVVW52  (7,5) 


are  functions  of  the  set  of  variables 


0^  02,  grad  01,  grad  02>  v,  grad  v 


(7.6) 


In  the  absence  of  the  stationary  rigid  solid,  invariance  requirements  under 
superposed  rigid  body  motions  would  indicate  that  the  velocity  v and  the  skew 
part  of  grad  v should  be  excluded  from  the  list  of  variables  (7.6).  However, 

r>J 

in  the  present  study  in  which  the  fluid  flows  through  a stationary  rigid  solid 
no  reduction  in  the  variables  (7*6)  is  possible. 

According  to  the  procedure  outlined  in  §3,  the  energy  equation  (7 • 3 ) is 
to  be  regarded  as  an  identity  for  all  motions  and  all  temperature  fields. 
Hence,  with  the  use  of  (7.4)  and  the  above  constitutive  assumptions,  applica- 
tion of  standard  techniques  results  in 


fi  = ^(0i,92)  + J"(0i,e2,v)  , 

^2  = ^(9l,e2)  + ^2(0l,e2,~)  » 

A 

* = ♦(e1,e2)  , p'ii  = + p2^2 


*i<Ve2»2.)  = 0 1 
9i»02’2^  = 0 * 

» p2^2  = 0 » 


p^l 


(7-7) 


(7.8) 


Hence,  end  depend  only  on  the  constituent  temperatures  end  coincide 

respectively  with  and  in  the  notation  of  §3-  Also,  with  the  help  of  (7.7) 
and  (7.8)  the  identity  (7.3)  reduces  to 
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A A 

(fc-  pi  aeT  «rod  e2  + p2  if  erad  V * ~ 


i 


A 


+ (£‘  Pl  17  Z)  • «rad  1“  Pi?!®!  ‘ P2?2®2 
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£ • grad.  61  - £ • grad  9g  - 0 


,,1 


(7-9) 


Next,  we  turn  to  the  inequality  (5.9)  with  given  by  (3.28)g  and  obtain 


A. 


^2 


Jt  - pi  3^sr8d  ®2-p2  ae^grad  ei  + E 


afj 


(7. 10  J| 


pi  *5r>«rad  0 » 


where  j*  and  T are  functions  of  the  variables  (7.6)  which  vanish  when  v and 
grad  v are  both  zero. 

In  the  rest  of  this  section,  we  limit  the  discussion  to  a linearized  theory 
about  an  equilibrium  state  in  which  v = 0 and  ®1=  9g  = 9,  9 being  a constant. 

We  also  assume  that  the  rigid  solid  is  homogeneous  and  isotropic  with  a center 
of  symmetry.  Thus,  we  set 


» 9(l+er+P)  , 92  ■ 9(l+a) 


(7.11)1 


and  retain  only  linear  terms  in  ar,g,v  in  all  constitutive  equations  for 


1 12 

£ >T,p  ,jp  However,  in  view  of  the  relations  (7.8)  and  the  energy 


identity  (7*9) > for  consistency  we  retain  both  the  linear  and  the  quadratic 
terms  in  a,p,v  in  the  expressions  for  , 5^ . Thus,  making  use  of 

(7.7)  and  remembering  that  the  rigid  solid  is  assumed  to  be  homogeneous  and 


isotropic,  we  write 


pA 


= Y0  - - Ygp  - i V3«  - yka?  - £ + Yv  • v , 


P2*2  - - 


v0  ‘ ^ia  • y2P  " * v3«2  - " Yv  • v , 


- „2 


(7.32)1 


# = 60  " 6la  " 52P  " * 63°2  “ “ * 65®2  ’ 


where  the  coefficients  6Q > • • - > 6^  are  related  to  yoi.,.ty  through  (7-7)y« 
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In  addition 


£ * + D2  grad  a + grad  g , 

p1  = d^v  - V grad  or  - Dg  grad  g , 

p2  = D?v  = Dg  grad  a - Dg  grad  g , 


(7.13) 


I f 
l £ 


T 

T = (jj*  + M-1)grad  v + (u  - U^Hgrad  v)  > 

Pi*i  "*Dio  ®“Dn  e+  pi*i  ’ 


(7.14) 


= D10  a f D11  ^ + p2^2  * 


where  y ,..*>V  Di»“*»Dll>  ^1  ar*  constarlts  whilc  5.1^2  ere  <5uadretic  func“ 


tions  in  the  variables 


a,  g,  grad  a,  grad  g>  v,  grad  ^ 


(7*15) 


and  the  partial  internal  production  of  entropies  %v%2  both  vanish  in  the 
equilibrium  state.  In  (7.l4)2>3  the  linear  parts  of  ?1,52  have  been  chosen 
so  that  the  identity  {(S)  is  satisfied  and  the  restrictions  on  the  forms  of 
7 or  f2  arise  through  the  quadratic  terms  in  (7*9) • Also,  by  (7.8),  (7‘H) 
and  (7.12),  the  expressions  for  the  partial  specific  entropies  are 

"A-  - i M * s2  + V*  65S  ’ 

0 

pA-'t  <1-4-  ' 

With  the  help  of  (/..«)]_  end  (Y.lh^,  the  inequality  (7.10)g  which  holds 
for  all  processes  may  be  used  to  show  that 

D=0  , D = 0 , D £ 0 , (j,SO  , - 0 1 

2 J 1 


(7.16) 


(7.17) 
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V-a.lt 


m 


Our  second  thermodynamic  condition  (5* 15 )> which  holds  in  equilibrium  with 
= rec*uces  the  present  problem  to 


X r 

- (£  + js " ) • grad  a > 0 when  0 - 0 , v » 0 
After  substitution  from  (7.13) » we  readily  obtain  the  condition 


(7.18) 


D5  + Dg  S 0 


(7.19) 


Further,  when  9^.  = ®2  0r  ^ = 0»  the  value  of  P*  “ p1«1  + p2«2  ia  60  + 61  + 63a  ae 
far  as  the  terms  linear  in  a are  concerned.  Then,  the  third  thermodynamic 
condition  (5 .19)  requires  that 


s3  a 0 


(7.20) 


Making  use  of  (7.1),  (7.13)n  (7.14)  and  (7.16)  and  retaining  only  the 

z o 


linear  terms,  the  entropy  balance  equations  become 


*4  |f  4 ‘5  at  - piBi  - V ■ D11S  + - 


(».'»4>  If  + It  0 PpBp4D,o“4Diie4V2“4VZ8  . 


(7.21)  1 
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whore  stands  for  the  Laplaclan  operator.  The  above  differential  equations, 
which  serve  to  determine  the  temperatures  a and  p,  are  independent  of  the 
velocity  v.  For  the  present  discussion,  we  do  not  record  the  differential 
equation  resulting  from  the  balance  of  momentum  which  is  an  equation  for  v 
and  p;  the  remaining  equation  involving  v is,  of  course,  the  incompressibility 


I 


condition  (7.1)^. 


Returning  to  (7. 21),  we  note  that  so  fR.r  we  have  only  obtained  the  thermo- 
dynamic restrictions  (7-19)  and  (7-20 ) on  the  constant  coefficients  which  occur  in 
the  entropy  balance  equations.  Our  previous  procedure  of  using  an  entropy  inequality 
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of  the  Clausius -Duhera  type  [see  the  inequality  (9.5)  in  the  paper  of  Craine,  Green 
and  Naghdi  (2)]would  yield  further  restrictions  on  the  coefficients  DC,D,,D0,IL, 
not  0n  neverti:iGleEIS>  we  not  adopt  this  approach  here. 

In  closing  this  section,  we  indicate  what  further  restrictions  may  be 
imposed  on  the  coefficients  6^>  ■ • • in  order  to  render,  under  suitable 
boundary  and  initial  conditions,  the  solutions  of  (7.21)  unique.  We  suppose 
that  the  stationary  rigid  solid  always  occupies  a region  of  space  R and  that 
the  fluid  also  occupies  the  same  region  so  that 


v . n = 0 on  dR  . (7.22) 

fsj 

1 o 

Recalling  (7. 13)  we  write  the  constitutive  equations  for  p ,p  in  the  forms 
2,3  **  ** 

p1  « D.  v + m1  , p2  = D7v  + m2  , (7.23) 

so  that 

11  2 2 

P • n - m • n , p^  • n = m • n on  dR  , (7.24) 

tsJ  r+J  sss  r**,  r*J 

in  view  of  (7*22). 

Consider  now  an  initial  boundary-value  problem  characterized  by  the  dif- 
ferential equations  (7.2l)  with  s =s  =0.  We  assume  that  the  dependent  variables 

2 2 1 

a,P  are:  (i)  of  class  C with  respect  to  x and  of  class  C with  respect  to  t in  the 
open  region  R for  tiO;  end  (ii)  of  class  C1  with  respect  to  x and  of  class  C with 
respect  to  t in  the  closed  region  R U dR  for  tSO.  In  addition,  we  assume  the 
Initial  and  boundary  data 
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RUdR  for  t £ 0 

a = 0 

. 3 

"■  0 

Oil 
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from  (7.21)  and  (7-25),  we  obtain 
34. 


(7.25) 


In  order  to  establish  uniqueness, 
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J ^(63a2  + 26UQfP+65P2)dv 

R 

= -J  ^D10a  + Dll^^  ' 2.1  • grad(a+p)  - m2  • grad  a]dv  , (7.26) 

ft 

which  is  an  identity.  Then,  with  suitable  classes  of  functions  a and  p,  the 
usual  type  of  energy  criteria  for  uniqueness  follow  from  (7*26)  provided 

63o2+26^+  6532  s.0  , 

(D10a  + Dll*)P  * 0 , (7.27) 

1 2 

-m  • grad(o+B)  - m • grad  a £ 0 , 

f*J  tst 

for  all  temperature  fields  a,ji.  Further,  with  the  help  of  (7.14)  and  (7.23), 
it  follows  from  (7.27)  that, 


VD8*°  , 06*o  , i*d6(d5  + d8)  s (d,i  + d6  + d9)2  , 
ho-0  • Dll*°  • »3*°  • 65SO  ’ ' 


(7.28) 


and  we  note  that  (7.28)^  and  (7*28)^  are  the  same  as  the  thermodynamic  restric- 
tions (7.19)  and  (7.20),  respectively. 

In  order  to  illustrate  some  of  the  features  of  the  equations  (7.21)  subject 
to  the  restrictions  (7,28),  we  consider  a simple  example  of  heat  conduction  in 
an  infinite  stationary  rigid  porous  solid  bounded  by  the  planes  x1  = 0,  x.^  * a 
and  containing  an  incompressible  viscous  fluid  at  rest.  At  time  t = 0,  the 
temperature  throughout  the  interior  of  the  solid  and  fluid  is  constant  and 
equl  to  0.  Let  the  temperature  of  the  walls  of  the  solid  be  changed  to,  and. 
Maintained  at,  another  constant  temperature.  Then,  when  0<x^<a,  the 
quantities  a,0  satisfy  (7.21)^  2 with  s^ = s2  - 0 and  are  functions  of  t and 
x^.  Also, 
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a - 0 , (3  = 0 (t'0  ; 0<x^<a)  , 


(7-29) 


a e a , 

o 


0 = 0 (x^  =>  0,a  ; t £ 0)  , 


where  aQ  is  a constant.  Following  a routine  method,  we  look  for  solutions  of 
(7.2l)]_  2 in  the  forms 

a ■ Ae^  sinCnnx^/e)  , 0 * Be^  sinCnTTx^/a)  , 

where  X satisfies  the  quadratic  equation 


(7.30) 


(6365  - 6^)X2  XC  C 63D6  + 6^(D5  + D8)  - 64(D5  + V D9  )}  9n  + 6^) 
+ (D6D8"D5D9)en+Dll(D5  + D8)0n  = 0 


0n  = n2rr2/a2  . 


Recalling  the  conditions  (7.28),  we  readily  obtain  the  restrictions 


63 65  - 62  S 0 , D*Da  - DqDQ  £ 0 , Dn n (Ds  + Dfl)  S 0 , 


J6V6  "59 


11 v "5  ' 8' 


63D6  + V°5  + “S'  ■ “l.(D5  + V V 8 0 • SDn  * 0 

and  then  we  conclude  that  the  roots  of  (7.33)  are  in  general  of  the  form 


X—On  * X="Pn  = » 


a > 0 , 0 > 0 , 

n ’ Kn 


where  «n>Pn  are  unequal  and  real,  or  are  of  the  form 


(7.31) 


(7-32) 


(7.33) 


(7.34) 


X = - orn  i 10^  (n-1,2,...)  , 


(7.35) 
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where  a, 8..  are  real  and  8 are  non-zero.  Exceptional  cases  can  arise  but 
n'm  rn 

these  can  be  dealt  with  similarly.  The  solution  to  the  initial-boundary 
problem  may  now  be  completed  in  a usual  way.  We  note,  in  particular,  one 
result  which  follows  from  the  solution,  namely  that 

tt-’OfQ  » B-*0  as  t - os  . 

According  to  (7.36),  which  is  an  expected  result,  the  tenperatures  in  the 
fluid  and  solid  tend  to  the  same  values  as  t -*  «*>  . 


(7.36) 
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8.  Further  remarks. 


We  close  this  paper  with  some  additional  remarks  concerning  the  present 
and  previous  procedures  of  studying  the  thermodynamics  of  mixtures  of  inter- 
acting continue  with  different  constituent  temperatures.  First,  we  observe 
that  the  field  equations  resulting  from  momenta  balances  for  each  constituent 
involve  only  the  mechanical  fields  while  those  obtained  from  balances  of  energy 
for  each  constituent  involve  both  the  thermal  and  the  mechanical  fields.  In 
the  traditional  developments  of  thermomechanics,  after  the  specification  of 
constitutive  equations,  the  former  field  equations  yield  a system  of  differential 
equations  for  the  determination  of  deformations  for  each  constituent  while  the 
latter  is  regarded  as  differential  equations  for  the  determination  of  the 
temperatures,  after  elimination  of  the  external  body  forces  with  the  use  of 
the  local  equations  of  motion.  In  contrast,  the  thermodynamical  development  of 
the  present  paper  is  such  that  the  balances  of  entropy,  which  involve  only  the 
thermal  fields,  provide  the  differential  equations  for  the  temperatures  in 
parallel  with  the  fact  that  the  differential  equations  for  the  deformations  are 
derived  from  the  balances  of  momenta.  Consistent  with  this,  the  energy  equation 
for  the  mixture  as  a whole  (after  the  elimination  of  the  external  body  forces 
and  the  external  heat  supplies)  is  employed  as  an  identity  for  all  thermo- 
mechanical processes;  and,  hence,  it  does  not  yield  any  new  differential  equa- 
tion for  temperatures  or  deformations.  It  is  noteworthy  that  the  energy  equa- 
tion for  the  mixture  as  a whole  as  used  here  provides  all  the  essential  results  ~ 
other  than  those  arising  from  inequalities  — concerning  the  relationships  between 
the  constitutive  response  functions  or  functionals  prior  to  any  appeal  to  a 
second  law  of  thermodynamics.  In  this  sense,  the  inequalities  (5*3)  or 
(5.8)  have  been  assigned  a more  subdued  role  relative  to  those  of  the  conservation  laws 

It  is  understood  that  the  differential  equation  for  temperatures,  as  ,ell  as 

those  for  the  deformations , contain  also  coupled  terms. 
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In  the  spirit  of  the  present  paper,  we  could  use  the  energy  equations  for 
each  constituent  as  energy  identities  provided  allowance  is  made  for  the  inter- 


action terms.  However,  this  can  he  accomplished  at  the  cost  of  introducing 
further  variables  requiring  constitutive  equations,  and  it  does  not  appear  likely 
to  yield  additional  useful  information. 

We  now  turn  to  an  entropy  inequality  of  the  Claus ius-Duhem  type  for 
mixtures,  namely 


v j » v « p r v r>  u ’ 11 

E dt  I PA  dV‘  Z J -T£dv+  S j da  £ 0 , 

am.  1 dt  a a ami  Jp  ea  a-1  9a 

a a a 

which  is  proposed  by  Truesdell  (ll)  , and  adopted  in  a number  of  recent  papers 
on  the  subject,  e.g.,  Green  and  Naghdi  (5,  6,  7)  and  Craine  et  al.  ( 2 ) . 

*•*-»  r*.*  /w 

The  inequality  (8.1)  may  be  regarded  as  a natural  generalization  of  the  Clauslus- 
Duhem  inequality  for  single  phase  continua  and  does  lead  to  reasonable  results 
in  some  special  situations.  For  regions  containing  no  surfaces  of  discontinuity, 
it  follows  from  (8.1),  (3.4)  and  (3*5)  that 
v v , D0^ 

A pA-J1^[eo<pa'Dt“tBaV-|>«r0talva  -E  •er*d  s0  <8‘2> 

for  all  thermo-mechanical  processes.  In  previous  papers  on  mixtures  the  external 
heat  supplies  were  eliminated  from  (8.2)  with  the  help  of  the  energy  equations 
postulated  for  each  constituent.  Such  a reduction  cannot  be  effected  here  since 
in  the  present  development  we  only  have  one  energy  equation.  As  a result, 
detailed  comparison  of  the  inequality  (8.2)  with  that  proposed  in  the  present 
paper,  namely  (5*9)>  is  difficult.  However,  some  comparison  can  be  made  when 
all  constituent  temperatures  9^  are  constrained  to  be  eqr  to  9(>0)  as  in 
(3.18).  Then,  using  the  energy  equation  (3*7)  and  the  various  definitions  in 
§3,  we  may  replace  (8.2)  by 


qa  • n 


(8.1) 
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+ E (u“  • va  + T®  • grad  va  - p®  . grad  0 ) 

. m /v  ^ 

0=1 


= E (p  w - pa  • grad  0 


a a ~ 


. v D«(VV) 
' " 2 [pJ Dt“ 


+ n (1)  -T)')}]  £ 0 . 

a 'a  'a 


(8.3) 


For  mixtures  whose  thermo-mechanical  response  is  such  that  ,T^  are 
independent  of  the  set  of  variables  (3.25),  the  specific  Helmholtz  free  energies 
and  the  specific  entropies  reduce  to  ^ and  then  (8.3)  becomes 


v v 

9 E p £ = E (p  w *p#  • grad  e)  £ 0 . (( 

o=l  * “ o=l  a a 

If  the  variables  * ,T1  ,T°,Lia,m  and  hence  p w in  (3-28).  do  not  depend  on  the 
1 a 'a  ~ o ro  a i 

temperature  gradient  grad  0 and  its  history,  then  it  readily  follows  from  (8.4) 

v 

that  E p w SO  and  this  is  the  same  as  our  inequality  (5-9)  • On  the  other 
o=l  a a v v 

hand,  provided  E p * grad  0<O,  the  inequality  (8.4)  could  allow  E p w to 
o=l  ~ o=l  a a 

be  negative  for  some  materials  undergoing  particular  admissible  processes. 


(8.4) 
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